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Abstract. We study a spatial branching model, where the underlying motion 
is Brownian motion and the branching is affected by a random collection of 
reproduction blocking sets called mild obstacles. We show that the quenched 
local growth rate is given by the branching rate in the 'free' region . When the 
underlying motion is an arbitrary diffusion process, we obtain a dichotomy for 
the local growth that is independent of the Poissonian intensity. Finally, and 
most importantly, we obtain the asymptotics (in probability) of the quenched 
(when d < 2) and the annealed (arbitrary d) global growth rates, and identify 
subexponential correction terms. 
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1. Introduction 

1.1. Model. The purpose of this note is to study a spatial branching model with 
the property that branching only takes place in a certain random region where it is 
spatially (and temporarily) constant. More specifically we will use a natural model 
for the random environment: mild Poissonian obstacles. 

Let w be a Poisson point process (PPP) on M"* with intensity v > and let P de- 
note the corresponding law. Furthermore, let a, /? > fixed. We define the branch- 
ing Brownian motion (BBM) with a mild Poissonian obstacle, or the '(i^, /?, a)-BBM' 
as follows. Let K denote the random set given by the a-neighborhood of uj: 

Then if is a mild obstacle configuration attached to uj. This means that given w, 
we define as the law of the (strictly dyadic) BBM on K'*, d > 1 with spatially 
dependent branching rate r{x,uj) :— (31k^{x). An equivalent (informal) definition 
is that as long as a particle is in K'^, it obeys the branching rule with rate /3, 
while in K its reproduction is completely suppressed and it does not branch at 
all. (We assume that the process starts with a single particle at the origin.) The 
process Z under is called a BBM with mild Poissonian obstacles. The total mass 
process will be denoted by \Z\. Further, W will denote d-dimensional Brownian 
motion with probabilities {P^;, x G W^}. Finally B{z,r) will denote the open ball 
centered at z g R'^ with radius r and [zj will denote the integer part of z G M: 
z := max{n G Z | n < z}. 

Remark 1 (Coupling). The following coupling between an ordinary "free" BBM 
and the one with mild obstacles may be useful to have in mind. Let us randomly 
"trim" the free historical BBM tree by looking at those branching points that are 
inside K and delete completely one of the two branches emanating from that point 
- the branch deleted is being chosen at random. This way we obtain a coupling: 
the randomly "trimmed" version of the free BBM has the same distribution as the 
BBM with mild obstacles. 

1.2. Motivation. Topics concerning the survival asymptotics for a single Brownian 
particle among Poissonian obstacles became a classic subject in the last twenty 
years, initiated by Sznitman, Bolthausen and others. It was originally motivated 
by the celebrated 'Wiener-sausage asymptotics' of Donsker and Varadhan in 1975 
and by mathematical physics, and the available literature today is huge — see the 
fundamental monograph jSz 98j and the references therein. 

In |E 2000| a model of a branching process in a random environment has been 
introduced. There, instead of mild obstacles, we introduced hard obstacles and 
instantaneous killing of the branching process once any particle hits the trap con- 
figuration K. 

In |EdH 2002j the model was further studied and a similar model with 'indi- 
vidual killing rule' was also suggested. Individual killing rule means that only the 
individual particle hitting K is eliminated and so the process dies when the last 
particle has been absorbed. The difference between that model and the model we 
consider in the present article is that now the killing mechanism is even 'milder' 
as it only inhibits temporarily the reproduction of the individual particle but does 
not eliminate the particle. 
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Working with a BBM in a Poissonian random environment in IE 20001 lEdH 2002) 
turned out to be quite challenging and resulted in some unexpected and intriguing 
results. 

An alternative view on our setting is as follows. Arguably, the model can be 
viewed as a catalytic BBM as well — the catalytic set is then K'^. Catalytic spatial 
branching (mostly for superprocesses though) has been the subject of vigorous 
research in the last twenty years initiated by Dawson, Fleischmann and others — 
see the survey papers |K 2000| and |DF 2002| and references therein. In those 
models the individual branching rates of particles moving in space depend on the 
amount of contact between the particle ('reactant') and a certain random medium 
called the catalyst. The random medium is usually assumed to be a 'thin' random 
set (that could even be just one point) or another superprocess. 

In more complicated models 'mutually' or even 'cyclically' catalytic branching 
is considered. (See again |DF 2002| .') 

Our model is simpler than most catalytic models as our catalytic/blocking areas 
are fixed, whereas in several catalytic models they are moving. On the other hand, 
while for catalytic settings studied so far results were mostly only qualitative we 
are aiming to get quite sharp quantitative result. 

For the discrete setting there is much less work available, one example is |KS 2003) . 
In that paper the authors study a branching particle system on Z'^, whose branching 
is catalyzed by another autonomous particle system on Z'^. There are two types of 
particles, the A-particles ('catalyst') and the B-particles ('reactant'). They move, 
branch and interact in the following way. Let Na{x,s) and Nb{x,s) denote the 
number of A- [resp. S-Jparticles at x G Z'' and at time s € [0, oo). (Here all 
Na{x,0) and Nb{x,0) with x g Z'* are independent Poisson variables with mean 
fiA [resp. /is].) Every ^-particle (B-particle) performs independently a continuous- 
time random walk with jump rate Da (Db)- In addition a i3-particle dies at rate 
5, and, when present at x at time s, it splits into two particles in the next ds time 
units with probability (3Na{x, s)ds + o{ds). Conditionally on the system of the A- 
particles, the jumps, deaths and splitting of the _B-particles are independent. The 
authors prove that for large P there exists a critical S separating local extinction 
regime from local survival regime. 

A further example of the discrete setting is |AB 200(1) . 

It appears that our proposed model of a BBM with "mild" obstacles has biolog- 
ical merit to it. First, one immediately has the following two biological interpreta- 
tions in mind: 

(1) Migration with unfertile areas (Population dynamics): Population 
moves in space and reproduces by binary splitting, except at randomly 
located reproduction-blocking areas. 

(2) Fecundity selection (Genetics): Reproduction and mutation takes place. 
Certain randomly distributed genetic types have low fitness: even though 
they can be obtained by mutation, they themselves are unable to reproduce, 
unless mutation transforms them to different genetic types. In genetics this 
phenomenon is called ^fecundity selection\ 

One question of interest is of course the (local and global) growth rate of the pop- 
ulation. Once one knows the global population size, the model can be rescaled 
(normalized) by the global population size, giving a population of unit mass (some- 
what similarly to the fixed size assumption in the Moran model or many other 
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models from theoretical biology) and then the question becomes the shape of the 
population. 

In the population dynamics setting this latter question concerns whether or not 
there is a preferred spatial location for the process to populate. In the genetic 
setting the question is about the existence of a certain kind of genetic type that is 
preferred in the long run that lowers the risk of low of fecundity caused by mutating 
into less fit genetics types. 

Of course, the genealogical structure is a very exciting problem to explore too. 
For example it seems quite possible that for large times the 'bulk' of the population 
consists of descendants of a single particle that decided to travel far enough (rcsp. 
to mutate many times) in order to be in a less hostile environment (resp. in high 
fitness genetic type area), where she and her descendants can reproduce freely. 

For example, a related phenomenon in marine systems (personal communication 
with Chris Cosncr) is when hypoxic patches form in estuaries because of stratifica- 
tion of the water. The patches affect different organisms in different ways but are 
detrimental to some of them. They appear and disappear in an effectively stochas- 
tic way. This is an actual system that has some features that correspond to the 
type of assumptions built into our model. 

It appears (personal communication with Bill Fagan) that a very relevant existing 
ecological context in which to place our model is the so-called "source-sink theory" . 
The basic idea is that some patches of habitat are good for a species (and growth 
rate is positive) whereas other patches are poor (and growth rate is zero or negative). 
Individuals c;an move between patches randomly or according to more detailed 
biological rules for behavior. 

Another kind of scenario where models such as the proposed one would make 
sense is in systems that are subject to periodic local disturbances (personal commu- 
nication with Chris Cosner). Those would include forests where trees sometimes fall 
creating gaps (which have various effects on different species but may harm some) 
or areas of grass or brush which are subject to occasional fires. Again, the effects 
may be mixed, but the burned areas can be expected to less suitable habitats for 
at least some organisms. 

Let us now return to our mathematical model and consider the following natural 

questions: 

(1) What can we say about the growth of the total population size? 

(2) What are the large deviations? (E.g., what is the probability of producing 
an atypically small population.) 

(3) What can we say about the local population growth? 

Of course, these questions make sense both in the annealed and the quenched sense. 

As far as the first question is concerned, recall that the total population of an 
ordinary (free) BBM grows in expectation (and almost surely) as e^^. To sec this 
note that in fact, for ordinary BBM, the spatial component plays no role, and hence 
the total mass is just a /3-rate pure birth process X. As is well known, the limit 
N limt^oo e~^*Xt exists a.s and in mean, and P(0 < N < oo) = 1. 

Turning to BBM with the Poissonian reproduction blocking mechanism, how 
much will the absence of branching in K slow the global reproduction down? Will 
it actually change the exponent /3? (We will see that even though the global re- 
production does slow down, the slowdown is captured by a sub-exponential factor, 
being different for the quenched and the annealed case.) 



BRANCHING BROWNIAN MOTION WITH "MILD" POISSONIAN OBSTACLES 



5 



Consider now the second question. Here is an argument to show the non-triviahty 
of the problem and to give a motivation. Let us ask the simplest question: what is 
the probability that there is no branching at all up to time t > 0? In order to avoid 
branching the first particle has to 'resist' the branching rate (3 inside K'^ . Therefore 
this question is quite similar to the survival asymptotics for a single Brownian 
motion among 'soft obstacles' — but of course in order to prevent branching the 
particle seeks for large islands covered by K rather then the usual 'clearings'. In 
other words, the particle now prefers to avoid the "Swiss cheese" K'^ instead of K . 
The second listed problem above is a possible generalization of this (modified) soft 
obstacle problem for a single particle. In fact, the presence of branching seems to 
bring genuinely new type of challenges into the analysis. 

Finally, the third question will be shown to be related to the recent paper 
|EK 2004] that treats branching diffusions on Euclidean domains. 

1.3. Expected global growth and dichotomy for local growth. Concerning 
the expected global growth rate we have the following result. 

Theorem 1 (Expected global growth rate). On a set of full V -measure, 



(1) S'"|Zt|=exp 



as t — > (X) 



(quenched asymptotics), and 

(2) (E®£:'^)|Zt| =exp[/3i-a(d,i^)t'*/('*+2)(l + o(l))], as i ^ oo 

( annealed asymptotics ), where 



c{d,v) Arf 



d + 2\ /2_A,^'^/^'^+') 
d 



andud is the volume of the d-dimensional unit ball, while Xd is the principal Dirich- 
let eigenvalue of —^A on it. 

Remark 2 (Interpretation of Theorem 1). Consider the annealed case first. Recall 
the coupling suggested in Remark ^ Having that coupling in mind, let us pretend 
for a moment that we are talking about an ordinary BBM. Then at time t one 
has e'^* particles with probability tending to one as t oo. For t fixed take a 
ball B = B{0,R) (here R = R{t)) and let K be so that B c K"^ (such a ball left 
empty by K is called a clearing). Consider the expected number of particles that 
are confined to B up to time t. These particles do not feel the blocking effect of K, 
while the other particles may have not been born due to it. 

Optimize R{t) with respect to the cost of having such a clearing and the prob- 
ability of confining a single Brownian motion to it. This is precisely the same 
optimization as for the classical Wiener-sausage. Hence one gets the expectation 
in the theorem as a lower estimate. 

One suspects that the main contribution in the expectation in ||2Jl is coming from 
the expectation on the event of having a clearing with optimal radius R{t). In other 
words, denoting by pt the probability that a single Brownian particle stays in the 
R{t)-hal\ up to time t, one argues heuristically that pte^* particles will stay inside 
the clearing up to time t 'for free' (i.e. with probability tending to one as t s- oo). 
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The intuitive reasoning is as follows. If we had independent particles instead of 
BBM, then, by a 'Law of Large Numbers type argument' (using Chebysev inequality 
and the fact that limt^ooPte'' = oo), roughly pte^* particles out of the total e^' 
would stay in the i?(i)-ball up to time t with probability tending to 1 as 1 1 oo. One 
suspects then that the lower estimate remains valid for the branching system too, 
because the particles are "not too much correlated" . This kind of argument (in 
the quenched case though) will be made precise in the proof of our main theorem 
by estimating certain covariances. 

The interpretation of the formula for the quenched case is different in the sense 
that large clearings (far away) are automatically (that is, P-a.s.) present; but it is 
similar in the sense that part of our job will be to show that inside such a clearing a 
large population is going to flourish (see the proof of Theorem 0] for more on this). 
O 

Concerning local population size we have the following (quenched) result. 

Theorem 2 (Quenched exponential growth). The following holds on a set of full 
P -measure: For any e > and any bounded open set % ^ B C. , 

( limsupe-('^-')*Zt(B) ^ oo \ > and ilimsnp e-^*Zt{B) < oo \ =1. 

Problem 3. What can one say about the distribution of the global and the local 
population size? Our theorems make it plausible that the global population size 
is a random multiple of exp [3t — c{d, v) Qpg*-)2/d (1 + o(l)) b& t go for almost 
every environment. This problem will be addressed in Section ll. 41 We will prove 
an appropriate formulation of the statement when the limit is meant in probability 
and d < 2. 

We now show how Theorem|21can be generalized for the case when the underlying 
motion is a diffusion. Let P be as before but replace the Brownian motion by an 
L-diffusion X on R"^, where L is a second order elliptic operator of the form 

— 1 2—1 

with aij,bi £ C"'{R"), a e [0,1), and the symmetric matrix aij{x) is positive 
definite for all a; G K*^. The branching i-diffusion with the Poissonian obstacles 
can be defined analogously to the case of BM. To present the result, we need an 
additional concept. Let 

Xc{L) = Xc{L, R'^) inf{A G K : 3u > satisfying {L - X)u = in M''} 

denote the generalized principal eigenvalue for L on M''. In fact Ac < 0, because 
il = 0. — see |P 1995| . section 4.4. 

The following theorem shows that the local behavior of the process exhibits 
a dichotomy. The crossover is given in terms of the local branching rate /3 and 
the term Xc{L): local extinction occurs when the branching rate inside the 'free 
region' K"^ is not sufficiently large to compensate the transience of the underlying 
L-diffusion; if it is strong enough, then local mass grows exponentially. Note an 
interesting feature of the result: the intensity v of the obstacles plays no role . 
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Theorem 3 (Quenched exponential growth/local extinction). Given the environ- 
ment uj, denote by the law of the branching L-diffusion. 

(i) Let 13 > — Ac(i) and let v > be arbitrary. Then the following holds on a 
set of full 'P -measure: For any e > and any bounded open set % ^ B d , 

P'^ ( limsupe("^"^<=(^)+")*Zt(B) =oo ) > 0. 
\ tToo / 

and 

( linisupe(-''-^'^(^»*Zt(B) < C50 ) =1. 

\ noo J 

(ii) Let (3 < —Xc{L) and let v > be arbitrary. Then the following holds on 
a set of full V -measure: For any bounded open set B C M'' there exists a 
P'^-a.s. finite random time to — tQ{B) such that Xt{B) = for all t > to, 
(local extinction). 

1.4. Quenched asymptotics of global growth. In this section we assume that 
d < 2 and investigate the behavior of the (quenched) global growth rate. Define 
the average growth rate by 

\og\Zt{Lo)\ 
rt = rt(uj) := . 

Replace now |Zt([x))| by its expectation Zf := E'^\Zt{uj)\, and define 

. . logZt 

rt = rt[uj) := — - — . 

Recall from Theorem ^ that on a set of full P-measure, 

(3) lim (log i)2/'^ (ft - (3) ^ ~c{d,iy). 

t — >oo 

We are going to show that an analogous statement holds for itself. 
Theorem 4. Let d < 2. On a set of full P-measure, 

(4) lim {\ogtf/'^{rt - (3)^ -c{d, v) in - probability. 

t — >oo 

(That is, loosely speaking, n ~ (3 — c{d, v)(\ogt)^'^/'^ .) 



1.5. Outline. The rest of this article is organized as follows. Section 2 discusses 
some further problems. Section 3 presents some preparatory results, while section 
4 contains the proofs of the theorems. Section 5 presents some additional problems 
that refer to certain parts of the proofs. In section 6 we discuss yet another result. 
Finally the appendix presents a proof of a statement that is needed in the proofs 
and that we consider to be of independent interest. 



2. Further problems 
In this section we suggest some further problems and directions for research. 
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2.1. More general branching. It should also be investigated, what happens 
when dyadic branching is replaced by a general one. In a more sophisticated popu- 
lation model, particles can also die — then the obstacles do not necessarily reduce 
the population size as they sometimes prevent death. 

For example critical branching requires an approach very different from the super- 
critical one, since taking expectation now does not provide a clue: E^\Zt{u))\ = 1, 
Mt > 0, Vw G Q. 

It is not even clear that 

P'^ (extinction) = 1 Vw e O. 

And if it is true, then what can we say about the asymptotic behavior (as t — > oo) 
of P(t > t)l (Here r denotes extinction time.) 

2.2. Super pro cesses with mild obstacles. A further goal is to generalize the 

setting by defining superprocesses with mild obstacles analogously to the BBM with 
mild obstacles. Recall the concept of an (L, /3, Q;)-superdiffusion: Let Mf denote 
the set of finite measures /x on M'' and let a, /3 denote functions in the Holder space 
satisfying a > and sup3,£][{ci l3{x) < oo. 

Notation 4. Let (X, , /i G A4[) denote the (L,f3,a) superdiffusion. That is, 
X is the unique Alf-valued (time-homogeneous) continuous Markov process which 
satisfies, for any bounded continuous : M'' i— > IR+ , 

(5) exp {Xt , -g) = exp ^- j ^ ij.{dx) u{t, x)j , 

where u is the minimal non-negative solution to 
d 

—u = Lu + (3u- aw^ on M** x (0, oo), 

(6) ^* 

(Here (f, /) denotes the integral J^^ i'{dx) f{x).) 

The definition of the superprocess with mild obstacles is straightforward: the 
parameter a vanishes on the set K and elsewhere it is positive (for example, a := 
uqIk' where ao is a positive constant). 

Similarly, one can consider the case when instead of a, the 'mass creation term' 
/3 is random, e.g. /3 := fiolx", Po > 0. Denote now by P'^ the law of this latter 
superprocess for given environment. We suspect that the superprocess with mild 
obstacles behaves similarly to the discrete branching process with mild obstacles 
when Xc{L + f3) > and P'^(-) is replaced by P'^(- | X survives). 

3. Preparations 
In this section we present some preparatory lemmas for the proofs. 
Lemma 5 (Expectation given by Brownian functional). Fix w. Then 



(7) E'^|Zt|=Eexp 



Jo 



ds 
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Proof. It is well known ('first moment formula' of spatial branching processes) 
that E^\Zt\ = (Ttl)(x), where u{x,t) :— {Ttl){x) is the minimal solution of the 
parabolic problem: 

du (\ 

In 

(8) u{-,Q) = 1, 

u > 0. 



Qa + /31^^ u on X (0, oo), 



This is equivalent (by the Feynman-Kac formula) with ■ 
We will also need the following result : 

Lemma 6 ( \EK 2004j . Theorem 3). Let Q < V , V he hounded from above 
and smooth. Let P denote the law of the {L^V) -branching diffusion X and let 
K ■.= \ciL + V). 

(i) Under P the process X exhibits local extinction if and only if Xc < 0. 

(ii) When Ac > 0, for any A < Ac and any open $ ^ B CC D, 

I limsupe-^*Xt(B) = oo I > and P^, i limsupe^^^* Xt{B) < oo ) =1. 

4. Proofs of Theorems 1-01 

4.1. Proof of Theorem 1. We can rewrite the equation ((T)) as 

E'^\Zt\= e/^'Ec^p - f P1k{Ws) ds. 

The expectation on the righthand side is precisely the survival probability among 
'soft obstacles', except that we do not sum the shape functions on the overlapping 
balls. However this does not make any difference with regard to the asymptotics 
(see |Sz 98| . Remark 4.2.2.). The statements thus follow from the well known 
Donsker-Varadhan type 'Wiener asymptotics' for soft obstacles ( jSz 98| Theorems 
4.5.1. and 4.5.3.; see also see |DV 75| ). ■ 

4.2. Proof of TheoremEJ Since in this case L = A/2 and since Ac(A/2, M"*) ^ 0, 
the statement immediately follows from Theorem |31 I 

4.3. Proof of Theorem Let g and G be two functions from C'''(M'') (here 
7 e (0,1]) with 

1b(0,o/2) < < ls(0,a) < G < 1, 

and let 

Fi(.) :=/3^5(.+a;), ^2(0 := ^ G(. + a>). 

Then 

(9) Vi < (31k^ < V2. 

Consider the operators L + Vi (i = 1 or i = 2) on M'* and let A^ = K:i^) denote 
their generalized principal eigenvalues. Note that Vi E C^{M.'^) and thus standard 
theory is apphcable. In particular, since Vi < V2, one has Aj < Ac (see Chapter 4 
in jP 1995p . As a first step, we show that in fact 

(10) Ac = Ac = Ac(L) + P for a.e. to. 
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By obvious comparison with the operator L + f3 on R'^, < Xc{L) + (3 for every 
ljj. On the other hand, one gets a lower estimate on X\ as foUows. Fix R > 0. 
By the homogeneity of the Poisson point process, for almost every environment, 
suppPi contains a clearing of radius R. Hence, by comparison, A;!. > X^^\ where 
A^^) is the principal Dirichlet eigenvalue of L + /3 on a ball of radius R. Since R can 
be chosen arbitrarily large and since limfl|oo A^^-* ~ Xc{L) + /3, we conclude that 
A^ > Ac(i) + /3 for almost every environment. 
From (jnj and ifTHI it follows that 

(11) A(L + /31kc) = Ac(L)+/3fora.e. w. 

The statements of the theorem follow from Hll|l along with Lemma El B 

4.4. Proof of Theorem |4l We give an upper and a lower estimate separately. 

Upper estimate: Let e > 0. Using the Markov and Jensen-inequalities along with 
the expectation formula QJ, we have that on a set of full P-measure: 



((logO'/'(n c{d, i^) > e) < e-'[(logt)^^^{E'^rt - (3) + c{d, ly)] 



o(l) as t — > oo. 



Lower estimate: We give a "bootstrap argument" : we first prove a weaker estimate 
which we will later use to prove the stronger result. The weaker result is as follows. 
Let {) < 6 < (3. Then on a set of full P-measure 

(12) lim P'^dZtl > e''*) = 1. 

t — >oo 

To prove (fT^ . recaU (fTHI and that Ac(A) = 0. Using these, take i? > large 
enough so that Ac ~ \c{Bji(Q)), the principal eigenvalue of the linearized operator 
+ (Hk" satisfies 

Xc{BRm > 5. 

Let be the process obtained from Z by introducing killing at dBfj{0). Then 



lim P{\Zt\ < e^*) < lim P(|Zf | < 



t — >OQ t — ^oo 



Let < <p — 4> he the Dirichlet eigenfunction corresponding to Ac on Bn{0), and 
normalize it by sup^g^^^g) 'Pi^) — 1- Then we can continue with 

< hm P((Zf,0) <e^*), 

where {Z^,(t)) := and {Z^''} is the Hth particle' in Z/^. Notice that 

Mt = Mf := e~^'^*(Z/^, (/)) is a non-negative martingale (see p. 84 in |EK 2004p . 
and define 

N := lim M*. 

t — >oc 

Since Xc{Bfi{0)) > S, the estimate is then continued as 

= lim P{Mt < e^-^"^* I iV = 0) P'^^\N = 0) < P^"'\N = 0). 

t — >oo 

We have that 

lim P{\Zt\ <e^') <P^"'\n = 0) 
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holds for all R large enough. Therefore, in order to prove ifT^ . we have to show 
that 

(13) lim p(^)(Ar > 0) = 1. 

i?,^oo 

Consider now the elliptic boundary value problem 

i Aw + (31k^ (u - ^ in Br{Q), 

^^^^ lima^^aSflXo) "(2^) = 0, 

u > in Bfl(O). 

The existence of a solution follows from the fact that Ac > by an analytical argu- 
ment given in jP 1996| . (Uniqueness follows by the semilinear maximum principle.) 

The argument below gives a probabilistic construction for the solution. Namely, 
we show that wr{x) := Px^\n > 0) solves (|Tl|l . To see this, let v = vr := 1 — wr. 
Let us fix an arbitrary time t > 0. Using the Markov and branching properties of 
Z at time t, it is straightforward to show that 



Since the left hand side of this equation defines a martingale in t, so does the right 
hand side. That is 

i 

defines a martingale. From this, it follows by Theorem 17 of |EK 2001] that v 
solves the equation obtained from the first equation of 1)1411 by switching u — v? to 
V? — u. Consequently, 'Wr{x) := Px^\n > 0) solves the first equation of ((T^ itself. 
That WR solves the second equation, follows easily from the continuity of Brownian 
motion. Finally its positivity (the third equation of l|14)) ') follows again from the 

fact that Ac > (see Lemma 6 in |EK 2004j 1. 

By the semilinear elliptic maximum principle (Proposition 7.1 in |EP99| : see also 
|P 1996p . wr{-) is monotone increasing in R. Using standard arguments, one can 
show that < w :— limj^^oo wr solves 1(11)1 too (see the proof of Theorem 1 in 

|p 

Applying the strong maximum principle to w := 1 — u;, it follows that w is either 
one everywhere or less than 1 everywhere. Now suppose that < w < 1. Then 

^Aw = /31i<-c(u;2 - ■a.) < 0. 

This contradicts the recurrence of the Brownian motion in one and two dimensions 
(see Chapter 4 in |P 19 95M . This contradiction proves that in fact w = 1 and 
consequently it proves ifT^ . 

We have now completed the first part of our "bootstrap" proof. 

Let us return to the proof of the lower estimate in Theorem 01 Let e > 0. We 
have to show that on a set of full P-measure, 

(15) lim P'^ ({logtf^'^irt - /3) + c{d, v) < -e) = 0. 

t—^OO \ J 
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To achieve this, we wiU define a particular function pt (the definition is given in 
(|24(l ) satisfying that as t ^ oo, 

t ( t 



(16) Pt = exp 



-c{d, v)- 



Using this function we are going to show a statement implying 115|) . namely, that 
for all e > there is a set of full P-measure, where 

(17) lim flog |Zf I <(3t + logpt - et{\ogt)-'^''^\ = 0. 

Let us first give an outline of the strategy of our proof. A key step will be in- 
troducing three different time scales, l{t), m{t) and t where i{t) = o{m(t,)) and 
m{t) = o{t) as t ^ oo. For the first, shortest time interval, we will use that there 
are "many" particles produced and they are not moving "too far away" , for the 
second (of length m{t) — t{t)) we will use that one particle moves into a clearing of 
a certain size at a certain distance, and in the third one (of length t — m{t)) we will 
use that there is a branching tree emanating from that particle so that a certain 
proportion of particles of that tree stay in the clearing with probability tending to 
one. 

To carry out this program, first recall the following fact (for example this can be 
found in the proof of Theorem 4.5.1 in |Sz 98\ ): Let 

l/d 



Arf / d 



c{d, y) V 



(recall that \d is the principal Dirichlet eigenvalue of — on the d-dimensional 
unit ball, and ujd is the volume of that ball) and let 

(18) p = p{i) := i?o(log^)^/'' - (loglog^)2, t > 0. 
Then, 

(19) P(3 4(w) > such that > £o(w) 3 clearing B{xo,p) with \xo\ <£) = !. 
Next, let £ and m be two functions M+ satisfying the following: 

(i) 

lim £(t) = oo, 



(ii) 

(iii) 
(iv) 

(v) 



lim 



t^oo \oge{t) 

£(t) ^ o{m{t)), as t - 
m(t) = o{£^{t)), as t - 



m{t) ^ o{t{\ogty^/'^), as i->oo. 

Note that (i) — (v) are in fact not independent, because (iv) follows from (ii) and 
(v). For example the following choices of £ and m satisfy (i) — (v): let £{t) and m{t) 
be arbitrarily defined for t G [0, e], and 

£(t) := ^i-i/Oogiogt)^ ^i-i/(2iogiogt)^ f^^, t>to>e. 
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Fix S e (0, (3) and define 

lit) := Lcxp(<5^(<))J. 
Let At denote the following event: 

At :={|Z,(,)|>/(t))}. 
By (I12|l we know that on a set of full P-measure, 

(20) lim P"(At) = 1. 

t — ^oo 

By (I20|l . for t fixed we can work on At C fl and consider I{t) many particles at 
time £{t). 

As a next step, we need some control on their spatial position. To achieve this, 
use Remark Q] to compare BBM's with and without obstacles, and then, use the 
following result taken from Proposition 2.3 in jEdH 2002] (the stronger, a.s. result 
is proved in Ky 2004| ). 

Let (only for the following statement) denote Z the BBM without obstacles 
starting at the origin with a single particle. Let R{t) — Usg[o,t] supp(Z(s)) denote 
the range of Z up to time t. Let 

(21) M{t) = inf{r > : R{t) C B{0, r)} for d > 1, 

be the radius of the minimal ball containing R{t). Then M{t)/t converges to \/2P 
in probability as t — > oo. 

Going back to the set of I{t) many particles at time £{t), H21() yields that even 
though they are at different locations, still for any e' > 0, with P"-probability 
tending to one, they are all inside the (\/2/? + e')i{t)-ha\l. 

Define 

p{t) := Ro[\og£it)]'/'' - [loglog^(t)]2, t > 0. 
Recall (|19|l . With P-probability one there is a clearing B = B{xQ,p{t)) such that 
1 2^0 1 < ^(i)j for all large enough t > 0. In the sequel we may assume without the 
loss of generality that |a;o| = £it). (Of course, xq depends on t, but this dependence 
is suppressed in our notation.) By the previous paragraph, with P'^-probability 
tending to one, the distance of xq from each of the I(t) many particles is at most 



{l + ^2l3 + e')£{t). 

Now, any such particle moves to B{xq, 1) in another m{t)—l{t) time with probability 
(7t, where (using {Hi) and {iv)) 

( [{l + V2B + e')t{t)f , ( [{l + ^+e')£(t)f \\ 
qt = exp — — — — ; h o — — — ^0 as t ^ oo. 

V 2[m{t)-£{t)] \ 2[m{t)~£{t)] J J 

Let the particle positions at time £(t) be zi, Z2, zj(^t) ^J^d consider the independent 
system of Brownian particles 

{W^,,; z = l,2,...,/(t)}, 

where ^^^,(0) = zf, i ^ 1, 2, /(t). In other words, {W^r, i = 1, 2, /(*)} just 
describes the evolution of the I{t) particles picked at time £{t) without respect to 
their possible further descendants and (using the Markov property) by resetting 
the clock at time £{t) . 

Let Ct denote the following event: 

Ct := {3i e {l,2,...,/(i)}, 30 < s < m{t) - £{t) such that e B{xoA)}- 
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By the independence of the particles, 
(22) hm sup P'^ {C^ I At) = hm sup(l 



)^(*) hmsup [(1 - qtf'^* 



qtHt) 



Since {Hi) implies that — o{i{t)) as i ^ oo and since (z) is assumed, one has 



(Z*e*^(*) = exp 



5l{t) + o{t{t)) 



oo as t 



2[m{t)^£{t)] 

In view of this, (|22|l imphes that linif^oo -P" (Cf I ^t) = 0. Using this along with 
(|20|l . it follows that on a set of full P-measure, 

(23) lim P"' (Ct) - 1. 

t — 'OO 

Once we know (|23f) . we proceed as follows. Recall that xq denotes the center of B 
and that P denotes the probability corresponding to a single Brownian particle W. 
Let cr^°'* denote the first exit time from B: 



ag°'' inf{s >0\Ws^B}. 



Abbreviate t* 
(24) 

where 



t — m{t) and define 



Pt 



sup 

xeB{xo,l) 



sup 

i;eS(0,l) 



C7B ■.^M{s>0\Ws^B{0,pit))}. 
It is easy to check that as i — > cxd, 

t* 



(25) 



Pt = exp 



-c(c?, i^) 



(l0g£(i))2/'i 

and using (m) and (w), it follows that in fact 

t 



t* 



(26) 



Pt = exp 



-c((i, i^) 



(logt)2/'i V(l0gt)2/'' 



A little later we will also need the following notation 

(27) pI:= sup P,(a^°'*>s)= sup P,(aB > s), 

xe-B(o,i) 



xeB(xo,l) 



With this notation, 



Pt 



P\^. 



By slightly changing the notation, let denote the BBM starting with a single 
particle at a; G i?; and let Z^'^ denote the BBM starting with a single particle at 
X ^ B and with absorbtion at dB (and still branching at the boundary at rate (3). 

Since branching does not depend on motion, \Z^'^ \ is a non-spatial Yule's process 
(and of course it does not depend on x) and thus for all x e B, 

(28) 



3N :^ lim e-^*|Zf-^| > 



almost surely (see Theorems III.7.1-2 in |AN 2004p . 

Recall Remark n By the coupling described there, it is clear that one can in 
fact define the process Z and the random variable N on the same probability space. 
Note that some particles of may re-enter B after exiting, whereas for Z^'^ that 
may not happen. Thus, by a simple coupling argument, one has that for all i > 0, 
the random variable |Zf (i3)| is stochastically larger than |Zf '^(i?)|. 



BRANCHING BROWNIAN MOTION WITH "MILD" POISSONIAN OBSTACLES 



15 



Recall that our goal is to show H17|l. and recall also (|24|l and (j26|l . In fact, we 
will prove the following, somewhat stronger version of H17|l : we will show that if 
the function 7 : [0, 00) [0, 00) satisfies linit^oo 7t = 0, then on a set of full 
P-measure, 



(29) lim (\Zt\ < -ft ■ e"' pt] ^ 0. 



't I ~^ It ' ^ 
Recalling t* ^ t — m{t), and setting 

7t := exp (^m{t) - e , for t > to > e, 

a simple computation shows that (|29|) yields H17|) . Note that this particular 7 
satisfies limt_>oo 7t = because of the condition (v) on the function m. 

By the comparison between \Zf{B)\ and \Z^'' {B)\ (discussed in the paragraph 
after (|28|l ) along with H23|) and the Markov property applied at time m{t), we have 
that 

lim P (\Zt\ < -ft ■ ef^^'pt) < lim sup P (\Z^:''iB)\ < jt ■ e^'' pt) . 

t-^QG \ / i— >00 \ / 

Consider now the J{x,t) := \Zp^\ many (correlated) Brownian paths starting 
a,t X ^ B and let us denote them by Wi, Wjj^, t). Let 

.J(x,t) 
i=l 

where 

A, := {W,{s) eB, V < s < 0- 
Then we have to show that 

(30) L := lim sup P (n^ < jt ■ e>^^' pt) = 0. 

Having the coupling between Z and N in mind, clearly, for all x ^ B, 



L = limt^oo sup^gs P ( < TT" j ^ 
(31) limt^oo sup.es [P < ipt) +P{N< 27t) 

Using the fact that lim^^oo 7t = and that N is almost surely positive, 

lim P{N < 27t) = 0; 

t — >oo 

hence it is enough to show that 

( \ \ 

(^2) /i"i!^^P(^(A^<2^*j=0- 

Let R denote the law of N and define the conditional laws 

P\-) :=P(. |7V = y), 2/>0. 

Then 

( ri^t 1 A f°° ^ , . „ f nf 1 



Define the conditional probabilities 
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where /it — ^t,y ■— l^^^^ J ■ Recall that H24|) defines pt by taking supremum over 



X and that \Zp^\ in fact does not depend on x. One has 



(33) 



P 



my) 



pv 



ye 



py 



<4^ 



As far as the second term of the integrand in ^6'6\ is concerned, the limit in H28|l 
implies that 



lim / R{Ay)py 



Zr|<^y = limP e 



-/3t* I yX,B I 



<-iV|^0. 



Let us now concentrate on the first term of the integrand in (|33ll . In fact, it is 
enough to prove that for each fixed K > 



(34) 



lim 

t—*oo 



1/K 



R{Ay) py 



ye 



I3t 



T<T:Pt]=0 



Indeed, once we know (jSU, we can write 



lim / R{dy)py 

t—>oo 



(35) 



lim 

t^OG 



Ridy) py 



l/K 



/ nf 1 

< 2^* 



R 



K 



R 



0, 



K 



Since this is true for all if > 0, thus letting K ] oo, 



lim / R{dy)py 



Returning to (|34|l . let us pick randomly /it many points out of the J{x,t) many 
particles — this is almost surely possible under py . (Again, 'randomly' means that 
the way we pick the particles is independent of their genealogy and their spatial 
position.) Let us denote the collection of these fit many particles by Mt- Define 

Then, one has 



(36) 



pv 



ye 



pt 



T<\pt]<P'' 



< 2^' 



We are going to use Chebysev's inequality and therefore we now calculate the 
variance. One has 

Var «) = fitipt - Pt) + Mf^t ~ l)^-^ '-^ , 

^J■t(^J■t - 1) 

where K{t,x) :— {{i,j) : i ^ J, 1 < *, j < /**}■ Now observe that 

^M^^l^f^^!^^ = iK^py)iAA,)-pl 

IJ'tyiJ't - tj 

where under P the pair is chosen randomly and uniformly over the fitifJ-t — 1) 
many possible pairs. 
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Let Q*'** and Q*-*^ denote the distribution of the death time of the most recent 
common ancestor of the ith and the jth particle under and under P, respectively. 
One has 

iI]®py){AA,)=pt [ [ p*_,>*)(0,s,dx)Q*'''(ds), 

Js=0 J B 

where 

p<*)(0,t,da;) := Po(W^t € da; | G S, z< t). 
By the Markov property applied at time s, 

P's / P^-s.xP^*HO,s,dx) =pt, 

and thus 

{E^py){A,A,)=pt f ^Q*'^(ds). 

Hence 

(37) Var'(«?) < - P^t) + f^tipt ~ l)p? • (/* - 1), 
where 

It := / b*]-iQ*'«(ds). 
Note that this estimate is uniform in x (see the definition oipt in (|24|l 'l. Define also 

/•CXD 

^t:- / [p*]-^Q(*)(ds). 

As a next step, we show that 

(38) lim Jt = 1. 

t — ^oo 

Since Jt > 1, thus it is enough to prove that 

lim sup Ji < 1 . 

t — ^CJO 

For r > we denote by A* := Xd^A, B{0,r)) the principal eigenvalue of on 
B{0, r). Since A* tends to zero as r t oo we can pick an i? > such that — A^ < /?. 
Let us fix this R for the rest of the proof. 

Let us also fix i > for a moment. From the probabilistic representation of 
the principal eigenvalue (see Chapter 4 in P 1995 ) we conclude the following: for 
e > fixed there exists a T(e) such that for s > T{i), 

logp* > (Ap(t) - e)s. 

Hence, for e > small enough (e < — A|j) and for all t satisfying Ap(t) > A|j + e 
(recall that limt^oo p{t) — oo) and s > T{e, t), 

(39) logp* >X*n-s. 

Note that T(e, t) can be chosen uniformly in t because^ limt^oo p{t) — oo, and so 
we will simply write T{e). Furthermore, clearly, T(e) can be chosen in such a way 
that 

(40) limT(e)=cx). 

no 

^In fact p{t) can be defined in a way that it is monotone increasing for large t's. 
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Depending on e let us break the integral into two parts: 



Jt=\ ^ [p*]-iQW(d5) + / [p^]-ig(')(d5) =: + jf 

Js=0 Js=T(f.) 



--T(e) 

We are going to control the two terms separately. 

Controlling J^^^ We show that 

(41) 3 lim J^^'^ < 1. 

t — ^oo 

First, it is easy to check that for all t > 0, (3'^*^(ds) is absolutely continuous, i.e. 
g(*)(ds) = 5(*)(s) ds with some 5^*) > 0. So 

l-Tie) ^ .Tie) 



=0 Js=0 



Evidently, one has [p*]^^ | as t ^ 00. Also, since Q''^\[a^b]) is monotone non- 
increasing in t for < a <h, therefore g^*\-) is also monotone non-increasing in t. 
Hence, by monotone convergence, 

lim = / 5(s)ds=lim / 5^*) (s)ds < lim / 5(*)(s)ds = l, 
Js=0 *^°°Js=o *-"°°Js=o 



t 

where g :— limt^oo 
Controlling Jj^^ Recall that 

(42) logp* > AJj • s, Vs > r(e). 

Thus, 

jf < / exp(-AJ^-5)Q(*)(d5). 

We will show that 
(43) 

lim lim / exp(-Ap • s) Q^*\ds) = lim lim / e-(^«+'')'' e'^'Q^'^ds) = 0. 

Recall that < /3 -I- A^. In order to verify (|43|) . we will show that given io > 
there exists some ^ < K — Kito) with the property that 

(44) g(*)(s) < Jfse^'^", for t > to, s e [to,*]- 

Indeed, it will then follow that 

lim lim / exp(-A|,-s)Q(*)(ds) = lim lim / exp(-A* • s) 5(*ns)(ds) 

eiO t^oo T(£)^oo t-»oo 



00 



<X lim / se~(^^+^)'' (ds) = 0. 

Recall that Q''*'* corresponds to the conditional law P(- | \Zf^\ > jit)- We now claim 
that we can work with P{- \ \Zf^\ > 2) instead of P(- | \ZP^\ > fit)- This is because 
if Qq^ corresponds to P(- | \Z^^\ > 2), then an easy computation reveals that for 
any e > there exists a to = io(e) such that for all t > to and for all < a < 6, 

g«([a, b]) ~ g(*'([a, b]) < 2(1 + e)Q^*)([a, b]); 
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thus, if 

(45) qI,'^ (ds) < Lse"^Ms on [to, t] for t > to 
holds with some L > 0, then also 

(46) Q(*' (ds) = 5^*' (s) ds < /sTse^'^Ms, t > to V £o, s € [^o, 

holds with K -.^ L + 2(1 + e). 

The bound (|45|l is verified in the appendix. 
It is now easy to finish the proof of l|38() . To make the dependence on e clear, let 
us write 4''^ = « = 1, 2. Then by one has that for all e > 0, 

(2) 

lim sup Jt < 1 + lini sup (e). 



Hence, (021 yields 



(2) 

lim sup Jt < 1 + lim lim sup (e) < 1, 

t— >oo eiO t^oc 



finishing the proof of H38() . 

Once we know (I38|l . we proceed as follows. Using Chebysev's inequality, one has 



By (|37ll . we can continue the estimate by 

r,2 



< 16 



Writing out fif , exploiting (|38|) . using that the lower limit in the integral is 1/K, 
and finally, dropping the —p^ term in the numerator, one obtains that 

(47) 

Ji/K \ye'^ 2 J 2 

(Recall that It in fact depends on y.) Since Yant^coPte^*' — oo, thus the first 
term on the righthand side of (|47|l tends to zero as t — > oo. Recall now that 

A*t •= L^^^T — J • ^ second term of (|47|l is concerned, it is easy to see that 

it also tends to zero as t ^ oo, provided one knows 

lim / R{dy){It-l) = Q. 



t — ^oo 



But Jjj°° R{dy){It ~1) — Jf — l and so we are finished by recalling (|38|l . Hence H34() 
follows. This completes the proof of the lower estimate in Theorem 0] ■ 



5. Some additional problems 



These problems were deferred to the section after the proofs, because the ques- 
tions themselves refer to parts of the proof. 
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Problem 7. The end of the proof for the lower esthnate in Theorem^lis basically 
a version of the Weak Law of Large Numbers. Using SLLN instead (and making 
some appropriate changes elsewhere), can one get 

liminf(logi)2/''(rt - P) > -c{d,r]) a.s. ? 

Problem 8. The question investigated in this paper was the (local and global) 
growth rate of the population. The next step can be the following: Once one knows 
the global population size \Zt\, the model can be rescaled (normalized) by \Zt\, 
giving a population of fixed weight. In other words, one considers the discrete 
probability measure valued process 

Then the question of the shape of the population for Z for large times is given by 
the limiting behavior of the random probability measures Zt, t>Q. (Of course, not 
only the particle mass has to be scaled, but also the spatial scales are interesting 
— see last paragraph.) 

Can one for example locate a unique dominant branch for almost every environ- 
ment, so that the total weight of its complement tends to (as t oo) zero? 

The motivation for this question comes from our proof of the lower estimate for 
Theorem 01 It seems conceivable that for large times the "bulk" of the population 
will live in a clearing within distance £(t) and with radius 

p{t) -.^ i?o[logf(t)]i/'^ - [log log ^(0]', t > 0, 

where 

lim £(t) = oo and lim = but lim ■; — ttt = 1- 

t^oo t-»oo t t— oo log lit) 

L. Mytnik asked how much the speed for free BBM reduces due to the presence 
of the mild obstacle configuration. As we have seen, for free BBM, the radius of the 
smallest ball covering the whole population grows linearly (the velocity is \/2f3). 
Since l{t) = o{t), it would be interesting to know if the speed in fact becomes 
sublinear. Note that in the above discussion about the shape we were only talking 
about the bulk of the population and not about individual particles travelling to 
very large distances from the origin. 

6. Annealed asymptotics of global growth 

In the annealed case we note that the following theorem can be justified by a 
method that is similar but much simpler than the one used to prove Theorem 0] 
Let Tt be as in Theorem 0] 

Theorem 5. For all d> 1, 

lim t^ {rt -13) = -c{d, v) in (E (g) P) - probability. 

t — >oo 

(That is, loosely speaking, rt ~ (3 — c{d, v) i^^/^'^+^^.j 

The upper estimate goes exactly the same way as in the proof of Theorem 0] by 
using what we know about the expectation. The lower estimate however becomes 
much easier. 

Now one does not need the two step 'bootstrap' argument as in the quenched 
case, and so there is no problem with higher dimensions (the d < 2 assumption 
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was crucial in the proof for the quenched case in the first part of the bootstrap 
method). The only remaining task is to show that the Law of Large Numbers is 
in force for the particle number inside the clearing around the origin. This can be 
done essentially the same way as in the last part of the proof of the quenched case. 
(The radius of the clearing is the same as in the classical annealed problem for a 
single Brownian particle, cf. the interpretation given for Theorem 

7. Appendix: Proof of the bound (|45|) 



We now give the proof of the bound (|45|l . In fact we prove a precise formula for 
the distribution of the death time of the most recent common ancestor, which, we 
believe, is of independent interest. The result and its proof are due to W. Angerer 
and A. Wakolbinger (personal communication). 

For simplicity we set (3—1; the general case is similar. Let us fix i > 0. Then 
for < M < i, one has 

(48) Qi'\s<u)^ (l-e-^)(Le--)^ ' 

and so the density is 

(,) _dQ* ^ _ ^ e-^ju -2+{u + 2)e-") + e-*(l - 2^e-" - e'^") 
/ W - W-^ (i_e-*)(l-e-«)3 

where I denotes Lebesgue measure on [0,t]. 

Proof of (|48|) : Consider the Yule population Yt :— \Z^^\ and recall that Qq'-* 
corresponds to P{- \ Yt > 2). The first observation concerns the Yule genealogy. 
Let us pick a pair of individuals from the Yule population at time t, assuming that 
Yt — j, j > 2. Denote by / the size of the population just before the coalescence 
time s of the two ancestral lines (where 'before' refers to backward time). That is, 
let / :— Fs+dt- Using some formulae from |EPW05j . we now show that 

(49) ^(^-^)-^tAv-^- 

J — 1 [1 — 1)1 I + 1 

Indeed, the distribution of / equals the conditional distribution of F given F < j, 
where F is defined as follows. First, the pure birth process K = (Ki) (with respect 
to 'Yule-time' i) is defined in Section 3.5 of (EPW05] . and setting n — 2, 

P{K,^, ^1\K, = 2) = 

i{i-l) 

(see formula (4.10) of the paper). Then the 'hitting time' F is defined (in the same 
section) by F := min{/ : Ki = 2}. The formula for the distribution of F (formula 
(2.3) of the paper) now becomes 

P{F <i) = P{K^ = 2) = i > 2. 

i + l 

Let i<j. Then 

P{I ^i)^P{F = i\ F < j) = P{K,-i = 1, X, = 2 I F < j) 

^P{K,^i = l,K, = 2\K,=2) 
^ P{K,^i = 1,K, = 2) 
P{K, = 2) 
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From the last three displayed formulae one arrives immediately at H49|l . 

Let us now embed the 'Yule time' into real time. Since a Yule population stem- 
ming from i ancestors has a negative binomial distribution, therefore, using the 
Markov property at times u and u + du, one can decompose 

(50) P{Yu = i-l, Yu+du ^i,Yt= j) = pi • p2 ■ Pa, 

where 

P2 — [i ^ 1) du and 

That is, 

P(Yu = i - l,Yu+du = t,Yt = j) 

= (i - 1) 2 e-(*-")Xl - e-(*-"')^"'e^"(l - e"")'-^ du. 

Since the pair we have chosen coalesce independently from the rest of the popula- 
tion, the random variables s and / are independent. Using the definition of s first 
and then the independence remarked in the previous sentence, and finally H49|l and 
(ISOl, 

P{se[u,u + du],Ys+dt = i,Yt^j) = P{l ^i,Yu = i-l,Yu+du = i,Yt= j) 

- P{I^i)P{Yu^i-l,Yu+du = t,Yt=j) 

i~2j + 
(1 - e-(*-"))^'-'e-"(l - e"")'-2 

for < u < 

Now, summing from j = i to cx), and from i = 2 to cx), and then dividing by 
P{Yt > 2) = 1 — e^*, one obtains (after doing some algebra) that for < u < t, 
r^Wt , J ^^ >^ _„2(2e-(*-")+i-l)(l-e-")'-2 

(1-e *)i(i + l) 

„ e-"(w-2 + (w-H2)e-")-he-*(l-2we-"-e~2«) 

= 2 • — -77:; 73 du . 

(1 - e-*)(l - e-«)3 

Equivalently, in integrated form, one has (|48|l . ■ 
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